3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression.

Enhancer hijacking, caused by structural alterations on chromosomes as well as extrachromosomal DNA (ecDNA), is a common cancer driver event. The complexity and ubiquity of structural alterations in cancer genomes make it difficult to identify enhancer hijacking using genome sequencing alone. Here we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking caused by structural alterations. HAPI analysis of HiChIP data from 34 cancer cell lines identified novel enhancer hijacking events that involve chromosomal rearrangements and activate both known and potentially novel oncogenes such as MYC, CCND1, ETV1, CRKL, and ID4, which we validated using CRISPRi assays and RNA-seq analysis. Furthermore, we found that ecDNAs often contain multiple oncogenes from different chromosomes, which causes nested enhancer hijacking among them. We found that ecDNAs containing MYC often harbor additional oncogenes from other chromosomes such as CDX2, ERBB2, or CD44 that co-opt MYC's enhancers for their overexpression, which we validated using dual-color DNA FISH and CRISPRi assays. These enhancer hijacking events involving multiple oncogenes on ecDNAs have important implications for therapeutic strategies that either target the co-opting oncogenes or the hijacked enhancers. Our publicly available HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation caused by chromosomal and extrachromosomal structural alterations.

Soluble factors in malignant ascites promote the metastatic adhesion of gastric adenocarcinoma cells.

BACKGROUND: Adenocarcinoma of the proximal stomach is the fastest rising malignancy in North America. It is commonly associated with peritoneal accumulation of malignant ascites (MA), a fluid containing cancer and inflammatory cells and soluble proteins. Peritoneal metastasis (PM) is the most common site of gastric cancer (GC) progression after curative-intent surgery and is the leading cause of death among GC patients. METHODS/RESULTS: Using a panel of gastric adenocarcinoma cell lines (human: MKN 45, SNU-5; murine: NCC-S1M), we demonstrate that prior incubation of GC cells with MA results in a significant (> 1.7-fold) increase in the number of cells capable of adhering to human peritoneal mesothelial cells (HPMC) (p < 0.05). We then corroborate these findings using an ex vivo PM model and show that MA also significantly enhances the ability of GC cells to adhere to strips of human peritoneum (p < 0.05). Using a multiplex ELISA, we identify MIF and VEGF as consistently elevated across MA samples from GC patients (p < 0.05). We demonstrate that agents that block the effects of MIF or VEGF abrogate the ability of MA to stimulate the adhesion of GC cells to adhere to human peritoneum and promote both ex vivo and in vivo metastases. CONCLUSION: Agents targeting MIF or VEGF may be relevant to the treatment or prevention of PM in GC patients.

To bind or not to bind: Cistromic reprogramming in prostate cancer.

The term "cistrome" refers to the genome-wide location of regulatory elements associated with transcription factor binding-sites. The cistrome of key regulatory factors in prostate cancer etiology are substantially reprogrammed and altered during prostatic transformation and disease progression. For instance, the cistrome of the androgen receptor (AR), a ligand-inducible transcription factor central in normal prostate epithelium biology, is directly impacted and substantially reprogrammed during malignant transformation. Accumulating evidence demonstrates that additional transcription factors that are frequently mutated, or aberrantly expressed in prostate cancer, such as the pioneer transcription factors Forkhead Box A1 (FOXA1), the homeobox protein HOXB13, and the GATA binding protein 2 (GATA2), and the ETS-related gene (ERG), and the MYC proto-oncogene, contribute to the reprogramming of the AR cistrome. In addition, recent findings have highlighted key roles for the SWI/SNF complex and the chromatin-modifying helicase CHD1 in remodeling the epigenome and altering the AR cistrome during disease progression. In this review, we will cover the role of cistromic reprogramming in prostate cancer initiation and progression. Specifically, we will discuss the impact of key prostate cancer regulators, as well as the role of epigenetic and chromatin regulators in relation to the AR cistrome and the transformation of normal prostate epithelium. Given the importance of chromatin-transcription factor dynamics in normal cellular differentiation and cancer, an in-depth assessment of the factors involved in producing these altered cistromes is of great relevance and provides insight into new therapeutic strategies for prostate cancer.

A predominant enhancer co-amplified with the SOX2 oncogene is necessary and sufficient for its expression in squamous cancer.

Amplification and overexpression of the SOX2 oncogene represent a hallmark of squamous cancers originating from diverse tissue types. Here, we find that squamous cancers selectively amplify a 3' noncoding region together with SOX2, which harbors squamous cancer-specific chromatin accessible regions. We identify a single enhancer e1 that predominantly drives SOX2 expression. Repression of e1 in SOX2-high cells causes collapse of the surrounding enhancers, remarkable reduction in SOX2 expression, and a global transcriptional change reminiscent of SOX2 knockout. The e1 enhancer is driven by a combination of transcription factors including SOX2 itself and the AP-1 complex, which facilitates recruitment of the co-activator BRD4. CRISPR-mediated activation of e1 in SOX2-low cells is sufficient to rebuild the e1-SOX2 loop and activate SOX2 expression. Our study shows that squamous cancers selectively amplify a predominant enhancer to drive SOX2 overexpression, uncovering functional links among enhancer activation, chromatin looping, and lineage-specific copy number amplifications of oncogenes.

Chromatin Looping Shapes KLF5-Dependent Transcriptional Programs in Human Epithelial Cancers.

Activation of transcription factors is a key driver event in cancer. We and others have recently reported that the Krüppel-like transcription factor KLF5 is activated in multiple epithelial cancer types including squamous cancer and gastrointestinal adenocarcinoma, yet the functional consequences and the underlying mechanisms of this activation remain largely unknown. Here we demonstrate that activation of KLF5 results in strongly selective KLF5 dependency for these cancer types. KLF5 bound lineage-specific regulatory elements and activated gene expression programs essential to cancer cells. HiChIP analysis revealed that multiple distal KLF5 binding events cluster and synergize to activate individual target genes. Immunoprecipitation-mass spectrometry assays showed that KLF5 interacts with other transcription factors such as TP63 and YAP1, as well as the CBP/EP300 acetyltransferase complex. Furthermore, KLF5 guided the CBP/EP300 complex to increase acetylation of H3K27, which in turn enhanced recruitment of the bromodomain protein BRD4 to chromatin. The 3D chromatin architecture aggregated KLF5-dependent BRD4 binding to activate polymerase II elongation at KLF5 target genes, which conferred a transcriptional vulnerability to proteolysis-targeting chimera-induced degradation of BRD4. Our study demonstrates that KLF5 plays an essential role in multiple epithelial cancers by activating cancer-related genes through 3D chromatin loops, providing an evidence-based rationale for targeting the KLF5 pathway. SIGNIFICANCE: An integrative 3D genomics methodology delineates mechanisms underlying the function of KLF5 in multiple epithelial cancers and suggests potential strategies to target cancers with aberrantly activated KLF5.

Author Correction: Human somatic cell mutagenesis creates genetically tractable sarcomas.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

C3D: a tool to predict 3D genomic interactions between cis-regulatory elements.

MOTIVATION: The 3D genome architecture influences the regulation of genes by facilitating chromatin interactions between distal cis-regulatory elements and gene promoters. We implement Cross Cell-type Correlation based on DNA accessibility (C3D), a customizable computational tool that predicts chromatin interactions using an unsupervised algorithm that utilizes correlations in chromatin measurements, such as DNaseI hypersensitivity signals. RESULTS: C3D accurately predicts 32.7%, 18.3% and 24.1% of interactions, validated by ChIA-PET assays, between promoters and distal regions that overlie DNaseI hypersensitive sites in K562, MCF-7 and GM12878 cells, respectively. AVAILABILITY AND IMPLEMENTATION: Source code is open-source and freely available on GitHub (https://github.com/LupienLabOrganization/C3D) under the GNU GPLv3 license. C3D is implemented in Bash and R; it runs on any platform with Bash (≥4.0), R (≥3.1.1) and BEDTools (≥2.19.0). It requires the following R packages: GenomicRanges, Sushi, data.table, preprocessCore and dynamicTreeCut. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Publisher Correction: ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters.

This corrects the article DOI: 10.1038/ncomms7186.

The Transcriptional Repressor Polycomb Group Factor 6, PCGF6, Negatively Regulates Dendritic Cell Activation and Promotes Quiescence.

Pro-inflammatory signals provided by the microenvironment are critical to activate dendritic cells (DCs), components of the innate immune system that shape both innate and adaptive immunity. However, to prevent inappropriate immune activation, mechanisms must be in place to restrain DC activation to ensure DCs are activated only once sufficient stimuli have been received. Here, we report that DC activation and immunogenicity are regulated by the transcriptional repressor Polycomb group factor 6 (PCGF6). Pcgf6 is rapidly downregulated upon stimulation, and this downregulation is necessary to permit full DC activation. Silencing PCGF6 expression enhanced both spontaneous and stimulated DC activation. We show that PCGF6 associates with the H3K4me3 demethylase JARID1c, and together, they negatively regulate H3K4me3 levels in DCs. Our results identify two key regulators, PCGF6 and JARID1c that temper DC activation and implicate active transcriptional silencing via histone demethylation as a previously unappreciated mechanism for regulating DC activation and quiescence.

Noncoding somatic and inherited single-nucleotide variants converge to promote ESR1 expression in breast cancer.

Sustained expression of the estrogen receptor-α (ESR1) drives two-thirds of breast cancer and defines the ESR1-positive subtype. ESR1 engages enhancers upon estrogen stimulation to establish an oncogenic expression program. Somatic copy number alterations involving the ESR1 gene occur in approximately 1% of ESR1-positive breast cancers, suggesting that other mechanisms underlie the persistent expression of ESR1. We report significant enrichment of somatic mutations within the set of regulatory elements (SRE) regulating ESR1 in 7% of ESR1-positive breast cancers. These mutations regulate ESR1 expression by modulating transcription factor binding to the DNA. The SRE includes a recurrently mutated enhancer whose activity is also affected by rs9383590, a functional inherited single-nucleotide variant (SNV) that accounts for several breast cancer risk-associated loci. Our work highlights the importance of considering the combinatorial activity of regulatory elements as a single unit to delineate the impact of noncoding genetic alterations on single genes in cancer.

Modulation of long noncoding RNAs by risk SNPs underlying genetic predispositions to prostate cancer.

Long noncoding RNAs (lncRNAs) represent an attractive class of candidates to mediate cancer risk. Through integrative analysis of the lncRNA transcriptome with genomic data and SNP data from prostate cancer genome-wide association studies (GWAS), we identified 45 candidate lncRNAs associated with risk to prostate cancer. We further evaluated the mechanism underlying the top hit, PCAT1, and found that a risk-associated variant at rs7463708 increases binding of ONECUT2, a novel androgen receptor (AR)-interacting transcription factor, at a distal enhancer that loops to the PCAT1 promoter, resulting in upregulation of PCAT1 upon prolonged androgen treatment. In addition, PCAT1 interacts with AR and LSD1 and is required for their recruitment to the enhancers of GNMT and DHCR24, two androgen late-response genes implicated in prostate cancer development and progression. PCAT1 promotes prostate cancer cell proliferation and tumor growth in vitro and in vivo. These findings suggest that modulating lncRNA expression is an important mechanism for risk-associated SNPs in promoting prostate transformation.

Pre-neoplastic epigenetic disruption of transcriptional enhancers in chronic inflammation.

Chronic periodontitis (CP) is a chronic inflammatory disease independently associated with higher incidence of oral cavity squamous cell carcinoma (OSCC). However, the molecular mechanism responsible for this increased incidence is unknown. Here we profiled the DNA methylome of CP patients and healthy controls and compared to a large set of OSCC samples from TCGA. We observed a significant overlap between the altered DNA methylation patterns in CP and in OSCC, suggesting an emergence of a pre-neoplastic epigenome in CP. Remarkably, the hypermethylated CpGs in CP were significantly enriched for enhancer elements. This aberrant enhancer methylation is functional and able to disrupt enhancer activity by preventing the binding of chromatin looping factors. This study provides new insights on the molecular mechanisms linking chronic inflammation and tumor predisposition, highlighting the role of epigenetic disruption of transcriptional enhancers.

Physical activity and genetic predisposition to obesity in a multiethnic longitudinal study.

Physical activity (PA) has been shown to reduce the impact of FTO variation and obesity genetic risk scores (GRS) on BMI. We examined this interaction using a quantitative measure of PA and two adiposity indexes in a longitudinal multi-ethnic study. We analyzed the impact of PA on the association between 14 obesity predisposing variants (analyzed independently and as a GRS) and baseline/follow-up obesity measures in the multi-ethnic prospective cohort EpiDREAM (17423 participants from six ethnic groups). PA was analyzed using basic (low-moderate-high) and quantitative measures (metabolic equivalents (METS)), while BMI and the body adiposity index (BAI) were used to measure obesity. Increased PA was associated with decreased BMI/BAI at baseline/follow-up. FTO rs1421085, CDKAL1 rs2206734, TNNl3K rs1514176, GIPR rs11671664 and the GRS were associated with obesity measures at baseline and/or follow-up. Risk alleles of three SNPs displayed nominal associations with increased (NTRK2 rs1211166, BDNF rs1401635) or decreased (NPC1 rs1805081) basic PA score independently of BMI/BAI. Both basic and quantitative PA measures attenuated the association between FTO rs1421085 risk allele and BMI/BAI at baseline and follow-up. Our results show that physical activity can blunt the genetic effect of FTO rs1421085 on adiposity by 36-75% in a longitudinal multi-ethnic cohort.

ABC: a tool to identify SNVs causing allele-specific transcription factor binding from ChIP-Seq experiments.

MOTIVATION: Detection of allelic imbalances in ChIP-Seq reads is a powerful approach to identify functional non-coding single nucleotide variants (SNVs), either polymorphisms or mutations, which modulate the affinity of transcription factors for chromatin. We present ABC, a computational tool that identifies allele-specific binding of transcription factors from aligned ChIP-Seq reads at heterozygous SNVs. ABC controls for potential false positives resulting from biases introduced by the use of short sequencing reads in ChIP-Seq and can efficiently process a large number of heterozygous SNVs. RESULTS: ABC successfully identifies previously characterized functional SNVs, such as the rs4784227 breast cancer risk associated SNP that modulates the affinity of FOXA1 for the chromatin. AVAILABILITY AND IMPLEMENTATION: The code is open-source under an Artistic-2.0 license and versioned on GitHub (https://github.com/mlupien/ABC/). ABC is written in PERL and can be run on any platform with both PERL (≥5.18.1) and R (≥3.1.1) installed. The script requires the PERL Statistics::R module. CONTACT: mlupien@uhnres.utoronto.ca SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters.

Chromatin interactions connect distal regulatory elements to target gene promoters guiding stimulus- and lineage-specific transcription. Few factors securing chromatin interactions have so far been identified. Here, by integrating chromatin interaction maps with the large collection of transcription factor-binding profiles provided by the ENCODE project, we demonstrate that the zinc-finger protein ZNF143 preferentially occupies anchors of chromatin interactions connecting promoters with distal regulatory elements. It binds directly to promoters and associates with lineage-specific chromatin interactions and gene expression. Silencing ZNF143 or modulating its DNA-binding affinity using single-nucleotide polymorphisms (SNPs) as a surrogate of site-directed mutagenesis reveals the sequence dependency of chromatin interactions at gene promoters. We also find that chromatin interactions alone do not regulate gene expression. Together, our results identify ZNF143 as a novel chromatin-looping factor that contributes to the architectural foundation of the genome by providing sequence specificity at promoters connected with distal regulatory elements.

BioID identifies novel c-MYC interacting partners in cultured cells and xenograft tumors.

The BioID proximity-based biotin labeling technique was recently developed for the characterization of protein-protein interaction networks [1]. To date, this method has been applied to a number of different polypeptides expressed in cultured cells. Here we report the adaptation of BioID to the identification of protein-protein interactions surrounding the c-MYC oncoprotein in human cells grown both under standard culture conditions and in mice as tumor xenografts. Notably, in vivo BioID yielded >100 high confidence MYC interacting proteins, including >30 known binding partners. Putative novel MYC interactors include components of the STAGA/KAT5 and SWI/SNF chromatin remodeling complexes, DNA repair and replication factors, general transcription and elongation factors, and transcriptional co-regulators such as the DNA helicase protein chromodomain 8 (CHD8). Providing additional confidence in these findings, ENCODE ChIP-seq datasets highlight significant coincident binding throughout the genome for the MYC interactors identified here, and we validate the previously unreported MYC-CHD8 interaction using both a yeast two hybrid analysis and the proximity-based ligation assay. In sum, we demonstrate that BioID can be utilized to identify bona fide interacting partners for a chromatin-associated protein in vivo. This technique will allow for a much improved understanding of protein-protein interactions in a previously inaccessible biological setting. BIOLOGICAL SIGNIFICANCE: The c-MYC (MYC) oncogene is a transcription factor that plays important roles in cancer initiation and progression. MYC expression is deregulated in more than 50% of human cancers, but the role of this protein in normal cell biology and tumor progression is still not well understood, in part because identifying MYC-interacting proteins has been technically challenging: MYC-containing chromatin-associated complexes are difficult to isolate using traditional affinity purification methods, and the MYC protein is exceptionally labile, with a half-life of only ~30 min. Developing a new strategy to gain insight into MYC-containing protein complexes would thus mark a key advance in cancer research. The recently described BioID proximity-based labeling technique represents a promising new complementary approach for the characterization of protein-protein interactions (PPIs) in cultured cells. Here we report that BioID can also be used to characterize protein-protein interactions for a chromatin-associated protein in tumor xenografts, and present a comprehensive, high confidence in vivo MYC interactome. This article is part of a Special Issue entitled: Protein dynamics in health and disease. Guest Editors: Pierre Thibault and Anne-Claude Gingras.

Enhancer alterations in cancer: a source for a cell identity crisis.

Enhancers are selectively utilized to orchestrate gene expression programs that first govern pluripotency and then proceed to highly specialized programs required for the process of cellular differentiation. Whereas gene-proximal promoters are typically active across numerous cell types, distal enhancer activation is cell-type-specific and central to cell fate determination, thereby accounting for cell identity. Recent studies have highlighted the diversity of enhancer usage, cataloguing millions of such elements in the human genome. The disruption of enhancer activity, through genetic or epigenetic alterations, can impact cell-type-specific functions, resulting in a wide range of pathologies. In cancer, these alterations can promote a 'cell identity crisis', in which enhancers associated with oncogenes and multipotentiality are activated, while those promoting cell fate commitment are inactivated. Overall, these alterations favor an undifferentiated cellular phenotype. Here, we review the current knowledge regarding the role of enhancers in normal cell function, and discuss how genetic and epigenetic changes in enhancer elements potentiate oncogenesis. In addition, we discuss how understanding the mechanisms regulating enhancer activity can inform therapeutic opportunities in cancer cells and highlight key challenges that remain in understanding enhancer biology as it relates to oncology.

Human somatic cell mutagenesis creates genetically tractable sarcomas.

Creating spontaneous yet genetically tractable human tumors from normal cells presents a fundamental challenge. Here we combined retroviral and transposon insertional mutagenesis to enable cancer gene discovery starting with human primary cells. We used lentiviruses to seed gain- and loss-of-function gene disruption elements, which were further deployed by Sleeping Beauty transposons throughout the genome of human bone explant mesenchymal cells. De novo tumors generated rapidly in this context were high-grade myxofibrosarcomas. Tumor insertion sites were enriched in recurrent somatic copy-number aberration regions from multiple cancer types and could be used to pinpoint new driver genes that sustain somatic alterations in patients. We identified HDLBP, which encodes the RNA-binding protein vigilin, as a candidate tumor suppressor deleted at 2q37.3 in greater than one out of ten tumors across multiple tissues of origin. Hybrid viral-transposon systems may accelerate the functional annotation of cancer genomes by enabling insertional mutagenesis screens in higher eukaryotes that are not amenable to germline transgenesis.

Laying a solid foundation for Manhattan--'setting the functional basis for the post-GWAS era'.

Genome-wide association studies (GWAS) have identified more than 8900 genetic variants, mainly single-nucleotide polymorphisms (SNPs), associated with hundreds of human traits and diseases, which define risk-associated loci. Variants that map to coding regions can affect protein sequence, translation rate, and alternative splicing, all of which influence protein function. However, the vast majority of sequence variants map to non-coding intergenic and intronic regions, and it has been much more challenging to assess the functional nature of these variants. Recent work annotating the non-coding regions of the genome has contributed to post-GWAS studies by facilitating the identification of the functional targets of risk-associated loci. Many non-coding genetic variants within risk-associated loci alter gene expression by modulating the activity of cis-regulatory elements. We review here these recent findings, discuss their implication for the post-GWAS era, and relate their importance to the interpretation of disease-associated mutations identified through whole-genome sequencing.

Variation at the DPP4 locus influences apolipoprotein B levels in South Asians and exhibits heterogeneity in Europeans related to BMI.

AIMS/HYPOTHESIS: Dyslipidaemia, a common feature of type 2 diabetes, is characterised by an increase in atherogenic particles, quantifiable through apolipoprotein B (ApoB) levels. Genetic studies of lipid levels have focused on Europeans; a study in South Asians could identify novel genes. METHODS: We tested 31,739 single nucleotide polymorphisms (SNPs) from ∼ 2,000 genes in 2,573 South Asians from the epidemiological arm of the Diabetes Reduction Assessment with Ramipril and Rosiglitazone Medication (DREAM) study (EpiDREAM) for association with ApoB and we tested two novel associations for replication in 1,181 South Asians from the INTERHEART case-control study. RESULTS: The SNP, rs4664443, within DPP4 was associated with ApoB (p = 7.98 × 10(-5)) in EpiDREAM. The observed association was replicated in the INTERHEART South Asians (one-sided p = 9.65 × 10(-3); combined two-sided p = 4.68 × 10(-6)). The rs4664443 SNP was not associated with ApoB among five other EpiDREAM ethnicities. However, because South Asians had a significantly lower mean BMI compared with other EpiDREAM ethnicities, we tested for and found an interaction between rs4664443 and BMI for ApoB among the Europeans, the largest subgroup in EpiDREAM (p = 4.14 × 10(-3) for interaction), observing an association with ApoB in Europeans with a BMI <25 kg/m(2) (p = 2.35 × 10(-3)), but not with a BMI ≥ 25 kg/m(2) (p = 0.21). The association between rs4664443 and ApoB among all EpiDREAM individuals with BMI <25 kg/m(2) was significant (n = 2,972; p = 1.44 × 10(-5)) compared with those with a BMI ≥ 25 kg/m(2) (n = 11,559; p = 0.81), and there was evidence of association among all genotyped individuals with a BMI <25 kg/m(2), including the INTERHEART South Asians (n = 3,601; p = 9.52 × 10(-7)). CONCLUSION/INTERPRETATION: Variation at the DPP4 locus is associated with ApoB in South Asians and displays heterogeneity related to BMI in other ethnicities.